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1. Introduction 

QCD is well established as theory of the strong interaction, and its perturbation theory expan- 
sion can be used to obtain quantitative predictions for observables in high energy particle colli- 
sions. QCD effects are omnipresent in hadronic collisions, and a detailed understanding of QCD is 
mandatory for the interpretation of collider data, for new physics searches and for precision studies. 
In this talk, I review the recent progress on applications of QCD at high energy colliders. 



2. Jets and Event Shapes 

Hadronic jets are the final state signatures of quark or gluon production in particle collisions 
at high energies. As such, they are important both as tool for precision studies of QCD, and in 
searches for new physics effects [1]. Jets are defined through a jet algorithm (a measurement and 
recombination prescription to reconstruct the jet momenta from measured individual hadron mo- 
menta). Jet algorithms must fulfill infrared-safety criteria, i.e. the reconstructed jet kinematics must 
be insensitive to radiation of soft or coUinear particles. Historically, two classes of jet algorithms 
were widely used at high energy colliders: cone-based algorithms and sequential recombination 
algorithms. Cone-based algorithms allow an intuitive understanding of the jets, and can be formu- 
lated in an infrared-safe manner [2]. Recombination algorithms are less intuitive, and their slow 
performance for a large number of final state particles was overcome only recently with the FastJet 
implementation [3]. Variants of these algorithms differ in the distance measure used to identify 
neighboring momenta, it turns out the the so-called anti-^r recombination algorithm results in per- 
fectly cone-shaped jets [4]. 

Much recent progress has been made recently in using jets as analysis tools. The concept of the 
jet catchment area [5] allows to obtain a geometrical interpretation of recombination algorithms, 
and to identify outside-jet regions, which can be used for underlying event studies. Aiming for 
the reconstruction of highly boosted massive particles, the study of jet substructure [6] has proven 
to be very promising. All decay products are first clustered in one fat jet, whose substructure is 
then resolved by lowering the resolution, resulting in a pronounced discontinuity once the particle 
decay is resolved. As one of the first results obtained using this procedure, the reconstruction of ttH 
(a reaction that could not be observed with standard cut-based methods due to the large standard 
model background) final states appears to become feasible [7]. Many more applications are under 
study. 

Closely related to jet observables are event shapes, which characterize the geometrical prop- 
erties of a hadronic final state. Distributions in several event shape variables were measured very 
extensively by LEP in view of precision studies of QCD. These results have a wide variety of appli- 
cations, ranging from precision measurements of ag, tests of resummation, study of hadronization 
effects, and tuning of multi-purpose Monte Carlo event generators. At hadron colliders, event 
shapes were only studied little up to now, and their definition is more involved due to the restricted 
final state region usually accessible in this environment. If defined properly, they can serve as 
tools for model-independent searches [8], and may be complementary to jet observables [9]. An 
extensive classification of event shapes at hadron colliders has been made very recently [8]. 
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3. Multiparticle Production at NLO 

The search for new physics signals at the CERN LHC will very often involve multi-particle 
final states, consisting of numerous jets, leptons, photons and missing energy. Quite in general, 
massive short-lived particles are detected through their decay signatures, as for example top quark 
pair production, which was first observed in final states with four jets, a lepton and missing energy. 

Meaningful searches for these signatures require not only a very good anticipation of the ex- 
pected signal, but also of all standard model backgrounds yielding identical final state signatures. 
Since leading-order calculations are affected by large uncertainties in their normalization and their 
kinematical dependence, it appears almost mandatory to include NLO corrections, which also al- 
low to quantify the jet algorithm dependence, and effects of extra radiation. For a long time, these 
corrections were available only for at most three final state particles. 

An NLO calculation of a w-particle observable consists of two contributions: the virtual one- 
loop correction to the w-particle production process, and the real radiation contribution from the 
{n + l)-particle production process. Both contributions are infrared divergent, and can be evalu- 
ated numerically only after extracting the infrared divergent contributions from the real radiation 
process. Several well-established and widely used methods exist for this task [10-17]. 

The evaluation of the one-loop multi-leg amplitudes poses a challenge in complexity (due 
to the large number of diagrams, and large number of different scales present) and stability (due 
to possible linear dependences among the external momenta). It has been known for long that any 
one-loop amplitude can be expressed as a linear combination of one-loop integrals with at most four 
external legs, plus a rational remainder. Enormous progress has been made in recent years in the 
systematic computation of the one-loop integral coefficients and rational terms. While previously 
established Feynman-diagram based techniques for tensor reduction and form factor decomposi- 
tion were successfully extended [18, 19] to multi-leg problems, a new arsenal of techniques was 
emerging from the use of unitarity and multi -particle cuts [20]. Using these, the one-loop integral 
coefficients of an amplitude can be inferred [21-24] without evaluation of all individual diagrams. 
An extension of these ideas is made by performing the reduction at the integrand level in the OPP 
method [25,26]. The rational coefficients can be determined in the same framework by extending 
the unitarity relations from four dimensions to higher-dimensional space-time [27-29]. 

Given the large number of different multi-particle final states of potential interest to new 
physics searches, an automation of NLO calculations is highly desirable. Based on existing multi- 
purpose leading order matrix element generators, the implementation of the real radiation contri- 
butions and their infrared subtraction terms is straightforward, and has been accomplished in the 
Sherpa [30], MadGraph [31-33] and Helac/Phegas [34] frameworks, as well as in the form of in- 
dependent libraries [35, 36], which complement already existing libraries in the MCFM [37, 38] 
and NLOJET-i~i- [39] packages. The automation of the virtual corrections is a much larger chal- 
lenge, which is currently being accomplished in several program packages based on the various 
available methods. A semi-numerical form factor decomposition is automated in the Golem pack- 
age [40]. Unitarity and multi-particle cuts are used in the BlackHat package [41], and the OPP 
method is automated in CutTools [42]. Numerical D-dimensional unitarity is applied in the Rocket 
package [43] and the Samurai package [44]; it also forms the basis of several currently ongoing 
implementations [45 , 46] . 
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Figure 1: Transverse momentum distribution of the third-hardest jet in Z" + 3] events at the Tevatron (left) 
and the LHC (right) to LO and NLO. Figure taken from [73]. 

Several NLO calculations of 2 — ;• 3 reactions at hadron colliders were completed recently. 
These include the production of two vector bosons and one jet [47-52], of a Higgs boson and two 
jets [53-56], of ttZ [57], and of three vector bosons [58-62]. Of a similar kinematical type are vec- 
tor boson fusion processes, which are computed to NLO accuracy in the VBFNLO package [63]. 
The current frontier of complexity are NLO calculations of 2 — )• 4 reactions. Several very important 
processes of this type have been computed most recently. 

An important channel for Higgs boson searches, and for subsequent determinations of Yukawa 
couplings, is the associated production of a Higgs with a heavy quark-antiquark pair, with the Higgs 
boson decaying into bb. The QCD background processes yielding ttbb final states were computed 
recently to NLO [64-67], displaying moderate but non-constant QCD corrections, which show a 
non-trivial dependence on the event selection cuts. A natural extension of these calculations are 
tt + 2j final states [68]. Extended Higgs sectors predict a sizable rate of associated production with 
bottom quark pairs, and the calculation of bbbb final states is in progress [69]. 

The final state signature of a vector boson and three hadronic jets is often relevant in generic 
new physics searches. NLO corrections of W + 3^ were obtained by two groups in the Rocket [70] 
and in the Blackhat-i-Sherpa [71,72] framework. The corrections to Z" + 3j were also obtained with 
Blackhat-i-Sherpa [73]. For both observables, corrections are moderate, and stabilize the QCD pre- 
diction to the ten per cent level, required for precision phenomenology, as can be seen in Figure |] 
from [73]. Knowledge of the NLO corrections to these processes allows many phenomenological 
studies, such as for example the stability of final state correlations [72] under perturbative correc- 
tions, and the optimal choice of scales in multi-scale processes [72-75]. A crossing of Z" + 3j is 
the process e^e^ — )• 5^, which was measured at LEP. The NLO calculation of it is in progress. 

4. Precision Observables at NNLO 



Few benchmark observables (e.g. jet cross sections, vector boson production, heavy quark 
production) are measured experimentally to an accuracy of one per cent or below. For a theoretical 
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interpretation of these observables, an NLO description (which has a typical residual uncertainty 
around ten per cent) is insufficient: extractions of fundamental parameters from these observables 
would be limited by the theory uncertainty. For a meaningful interpretation of these observables, 
NNLO corrections are mandatory. Likewise, NNLO corrections are required for a reliable descrip- 
tion of observables with potentially large perturbative corrections, like Higgs boson production. 

The calculation of NNLO corrections to an n-particle final state requires three ingredients: the 
two-loop matrix elements for the «-particle production, the one-loop matrix elements for the {n + 
1) -particle production and the tree-level matrix elements for (« + 2) -particle production. The latter 
two contributions develop infrared singularities if one or two particles become soft or collinear, 
requiring a subtraction method to extract these infrared poles, which are then combined with the 
virtual corrections to yield a finite prediction. The two major challenges of NNLO calculations are 
the two-loop matrix elements and the handling of the real radiation at NNLO. Up to now, two types 
of approaches to real radiation have been applied in NNLO calculations of exclusive observables. 
The sector decomposition method [76-78] is based on a systematic expansion in distributions, fol- 
lowed by numerical integration over many different small phase space sectors. Subtraction meth- 
ods search to approximate the full real radiation contribution by subtraction terms in all unresolved 
limits; these terms are then integrated analytically. While many subtraction methods have been 
worked out at NLO, only two methods have so far yielded results at NNLO: the antenna subtrac- 
tion method [79] for processes in e+e^ annihilation, and the ^j- -subtraction [80] for hadron collider 
processes in specific kinematic configurations. Alternative approaches are under intensive develop- 
ment [81-83]. A combination of subtraction with sector decomposition [84] may hold the potential 
to become a general multi-purpose method. 

The dominant Higgs boson production process is gluon fusion, mediated through a top quark 
loop. This process has been computed (in the infinite top mass limit) to NNLO accuracy in a 
fully exclusive form including the Higgs boson decay, i.e. allowing for arbitrary infrared-safe final 
state cuts, both using sector decomposition [85-87] and using ^j- -subtraction [88,89]. These re- 
sults can be directly applied to the Higgs boson search at the Tevatron, based on a neural network 
combination of many different kinematical distributions [90]. Finite top mass effects at NNLO 
were derived most recently [91-93] for the inclusive gluon fusion cross section. At this level of 
precision, mixed QCD and electroweak corrections [94] become equally important. The gluon fu- 
sion reaction can be mediated through loops involving any type of massive color-charged particles, 
thereby offering an indirect constraint on physics beyond the standard model, such as supersym- 
metric particles [95-99], extra heavy quark families [100] or color-octet scalars [101]. 

Another very promising Higgs discovery channel is vector boson fusion. The factorizable 
NNLO corrections to the inclusive cross section for this process are closely related to inclusive 
deep inelastic scattering. They were computed very recently [102], and turn out to be rather small, 
resulting in a high theoretical stability of the prediction. This channel can be equally sensitive on 
supersymmetric contributions [103]. 

Fully exclusive NNLO corrections to vector boson production have equally been derived using 
sector decomposition [104, 105] and with ^j- -subtraction [106], including the leptonic vector boson 
decay. Observables derived from vector boson production are very important for precision studies 
of the electroweak interaction, and for the determination of the quark distributions in the proton. 
Using the newly obtained results, the NNLO corrections (and their uncertainty) to the lepton charge 
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Figure 2: Lepton charge asymmetry at the Tevatron at NLO and NNLO, compared to CDF data. Figure 
taken from [107], 



asymmetry [107] can be quantified, see Figure §, and this observable can be consistently included 
into NNLO fits of parton distributions. 

Jet production observables have been computed to NNLO only for e^e^ annihilation up to 
now. Two implementations of the NNLO corrections to e^e^ — ;■ Sj and related observables are 
available [108-115], both based on antenna subtraction. The magnitude of the NNLO correc- 
tions differs substantially between different event shape observables; including these new NNLO 
corrections, LEP data on event shapes and jet cross sections were reanalyzed in view of an im- 
proved determination of the strong coupling constant. In general, an improved consistency among 
different observables was observed. To use measurements over an extended kinematical range, 
resummation of large logarithmic corrections in the two-jet limit is needed. This is available to 
next-to-leading logarithmic accuracy (NLLA) for all shape variables [1 16, 1 17], and to N^'LLA for 
thrust [118, 119] and heavy jet mass [120] distributions. The by-now limiting factor in precision 
physics with event shape observables in e^e^ annihilation is the description of the parton-to-hadron 
transition (hadronization), which was previously modeled from parton shower based event genera- 
tors. Substantial differences are observed between the LEP-era programs and more modern genera- 
tors, and to analytic approaches to hadronization, based on the shape function formalism [118-120] 
and on a dispersive model [121-123]. The recent determinations of the strong coupling constant 
from event shapes and jet cross sections at NNLO [118-120, 123-127] are summarized in Figure ||. 
Electroweak NLO corrections to jet observables [128-130] are potentially of the same numerical 
importance as NNLO QCD corrections, and could be included in future studies. 

In view of the very precise jet production data from HERA [131, 132] and the Tevatron [133, 
134], the derivation of NNLO corrections to jet cross sections in hadronic collisions is of high 
priority. The relevant two-loop matrix elements for hadronic collisions and for deep inelastic scat- 
tering [135] are known for some time already, and substantial progress is being made to extend 
the antenna subtraction method to include hadrons in the initial state. The proper functioning of 
this method on the gg — )• Ag subprocess to hadronic dijet production has been demonstrated [136] 
most recently. The integrated forms of all antenna functions have been derived for one parton in 
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Figure 3: Determinations of a^ from event shapes and jet cross sections in e+e^ annihilation at NNLO, 
compared to the Particle Data Group world average. Experimental errors are indicated in black, theoretical 
errors in blue. 



the initial state [137], the case of two initial state partons [138] is work in progress. 

The large number of top quark pairs expected to be produced at the LHC will allow for pre- 
cision top quark studies, requiring NNLO accuracy on the theoretical side. The relevant two-loop 
matrix elements were first derived in the high energy limit [139, 140]. The exact qq — >■ ff matrix ele- 
ment is known numerically [141], substantial parts of it have been confirmed by an analytic calcula- 
tion [142,143]. The one-loop self-interference contributions are equally known [144-146] The ma- 
trix elements with one and two extra partons form part of the tt + j production at NLO [147-149]. 
Methods to handle real radiation at NNLO in the presence of massive top quarks are currently un- 
der intensive development. Generalizing the subtraction method of [10] to NNLO and numerically 
integrating the relevant subtraction terms using sector decomposition [84] may provide a powerful 
method by combining the virtues of both approaches. 



5. Parton Distributions 



The parton distribution functions in the proton are a crucial ingredient in all hadron collider 
cross sections. They are determined (at LO, NLO or NNLO) from global fits [150] to a variety 
of data sets from fixed target experiments, from HERA and from the Tevatron. On the theory 
side, the parton distribution fits require the DGLAP splitting functions (which govern the scale 
evolution of the parton distributions, and are known to NNLO [151, 152]) and coefficient functions 
for each process considered in the global fit. These coefficient functions are known to NNLO for 
inclusive DIS [153], the Drell-Yan process [105] and for heavy quark production in DIS [154], but 
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only to NLO for jet production observables. The fit procedure must incorporate experimental and 
theoretical errors in a consistent manner. 

Global fits of parton distributions are performed by various collaborations, with slight dif- 
ferences in the methodology [150]. Recent sets of parton distributions are from MSTW [155], 
CTEQ [156], JR [157], NNPDF [158] and ABKM [159]. A comparison of them shows that the 
quark distributions are known rather precisely at large x, while the gluon distribution is uncertain 
to within ten per cent at large x, and systematic differences between the fits exist within errors. For 
small values of x < 10^^, uncertainties on the distributions become very large. 

6. Infrared Structure and Resummation 

The perturbative expansion of QCD observables in the strong coupling constant is reliable if 
only a single hard scale is present, it becomes problematic for observables depending on several 
hard scales, leading to large logarithmic corrections at all orders. In these cases, a rearrangement 
of the perturbative series by means of a resummation of large logarithmic corrections often appears 
more suitable. 

Resummation of leading logarithmic corrections is accomplished by event generators [160] 
based on parton showers, initially based on leading order calculations. Parton showers can be com- 
bined with fixed order NLO calculations in the MC@NLO [161] or the POWHEG [162] approach. 
The MC@NLO event generator already covers a large number of different processes, with W^t 
production [163] and H^^t production [164] among the most recent additions. Within POWHEG, 
single top production [165] and Higgs production in vector boson fusion [166] were accomplished 
most recently. The POWHEG box [167] provides users with a framework for implementing exist- 
ing NLO calculations in this framework. 

A detailed understanding of the infrared structure of QCD can be gained from the observation 
that infrared poles in loop amplitudes translate into large logarithms in real radiation processes 
and vice versa. This relation can be applied successfully in both directions: for example to pre- 
dict infrared poles at two loops from resummation [168, 169] and to extract large-x resummation 
constants [170] from the poles of the QCD form factors. By relating the infrared poles in QCD to 
ultraviolet poles in soft-coUinear effective theory (SCET) [171], it becomes possible to express the 
infrared pole structure of QCD amplitudes by a multiplicative renormalization in SCET. Based on 
constraints [172, 173] and symmetry arguments, it becomes possible to conjecture that the infrared 
pole structure of massless QCD multi-loop amplitudes is uniquely determined [172-175] by the 
cusp anomalous dimension and the coUinear anomalous dimensions of the external particles. 

The resummed description of an observable consists [176, 177] of a hard coefficient, a soft 
function, jet functions containing final state coUinear radiation and parton distributions containing 
initial state coUinear radiation. In SCET [171], each of these elements is identified with an operator 
or a non-local function. The resummation [178, 179] then proceeds by computing their anomalous 
dimensions and solving the respective evolution equations. First applications of SCET-based re- 
summation are the thrust [118, 119] and heavy jet mass [120] distributions in e+e^ annihilation, 
the inclusive Drell-Yan and Higgs production [178, 180] and inclusive photon production [181]. 
This topic is currently under fast development, any many yet open issues, like jet production and 
radiation off incoming partons [182, 183] are being addressed. 
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4: Top quark invariant mass distribution at the Tevatron in fixed order expansion (left) and resum- 
(right). Figure taken from [190]. 



Many of the constraints used to obtain the all-order conjecture for massless QCD amplitudes 
do not apply in the presence of particle masses. Consequently, the pole structure of massive am- 
plitudes is more involved; in particular, it contains multi-particle correlations [184], which were 
absent in the massless case. Only recently, a prediction of the infrared poles to two-loop order has 
been accomplished [184,185]. With these results, the resummation of the top quark pair production 
cross section to third logarithmic order (NNLL) could be completed. While dominant contributions 
at this order were known for some time [186, 187] the full corrections have been obtained now in 
two approaches: based massive soft anomalous dimensions [188, 189] and by using SCET [190]. 
The tt invariant mass distribution is compared in fixed order and resummed expansion in Figure 0, 
taken from [190]. It can be seen that the resummation has only moderate numerical impact on the 
central value, but results in a substantial reduction of the scale uncertainty. By expanding the re- 
summed results to fixed order, one can in turn approximate the NNLO corrections to the top quark 
production cross section [186, 191, 192]. 



7. IMulti-loop Results 

Several reasons motivate the derivation of perturbative corrections even beyond NNLO. Such 
theoretical accuracy is required to describe very precisely measured quantities (like sum rules or 
total decay rates) or for processes with a very slowly converging perturbative series. Moreover, 
these corrections allow insight into the infrared structure of QCD at high orders, and determine 
resummation coefficients. 

The current frontier of complexity for QCD loop amplitudes are 2 — )• 4 processes at one loop, 
2 — )• 2 processes at two loops, 1 — )• 2 processes at three loops and 1 — )• 1 processes at four loops. 
Various innovative techniques have allowed substantial progress on multi-loop QCD calculations 
in recent times: a substantial reduction of complexity is achieved by reducing the large number 
of integrals appearing in a calculation to a small number of master integrals by exploiting linear 
relations among the integrals [193]. Various techniques have proven successful in the derivation 
of these master integrals: for example the Mellin-Barnes transformation [194, 195], differential 
equations [196-199] in masses and momenta, and the sector decomposition technique [76,200]. 
The reduction to master integrals is usually based on a lexicographic ordering (the Laporta algo- 
rithm [201]). Implementations of this algorithm are available in several computer algebra frame- 
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works: the AIR package [202] in Maple, the FIRE package [203] in Mathematica, and the Reduze 
package [204] as a stand-alone C++ implementation based on Ginac and CLN. 

The Mellin-Barnes method allows to express master integrals in a systematic manner [205] in 
a form suitable for analytical or numerical evaluation. Recent results obtained with this method 
include the massless three-loop QCD form factor integrals [206-208]. Another very successful 
method is the glue-and-cut technique, which exploits relations among topologically different mas- 
ter integrals. It recently led to the derivation of the four-loop propagator master integrals [209]. 
Many of these results were validated independently [210] using the sector decomposition tech- 
nique. 

The quark and gluon form factors are the simplest infrared-divergent objects in QCD. They 
allow the determination of resummation coefficients and enter benchmark reactions like Drell- 
Yan and Higgs production. They were recently computed to three loops [211, 212]. The static 
QCD potential formed by an infinitely heavy quark-antiquark pair is an important ingredient in 
the determination of heavy quark masses from sum rules, and in the description of top quark pair 
production at threshold. It was computed at three-loop accuracy most recently [213,214]. 

Many important observables can be expressed as two-point functions: total decay rates, sum 
rules and moments of structure functions. Massless two-point functions were obtained recently 
at four loops. Among the results derived in this context are the hadronic /?-ratio and the T-decay 
rate [215], which allow some of the most accurate determinations of the strong coupling constant: 
a,{Mz)R = O.II90± 0.0026 and as{Mz)x = 0.1202 ±0.0019. Most recent results [216] are the 
polarized Bjorken sum rule, the Adler function and the Crewther relation. 

8. Conclusions 

QCD is crucial for the success of the LHC physics programme in understanding signals and 
backgrounds, knowing parton distribution functions, and using jets and event shapes as analysis 
tools. Particle theory is getting ready for this challenge on many frontiers: with improved jet 
algorithms and event shape definitions, with an enormous progress on NLO calculations for multi- 
leg final states, with first NNLO results for precision observables, with an emerging understanding 
of the all-order structure of infrared singularities, and with landmark results at three loops and 
beyond. 
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